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ABSTRACT 

We have mapped the 2.6 mm CO J = 1 ^ emission in three optically selected "Toomre Sequence" 
mergers (NGC 520, NGC 3921, NGC 4676). The molecular gas distribution is well resolved by the 
observations. For NGC 520 and NGC 4676A, the nuclear gas concentrations form a disk- or a ring-like 
structure, and the gas kinematics are regular and are consistent with simple rotation. Discrete molecular 
gas complexes are found along the stellar bar in NGC 4676B, and the gas kinematics is consistent with 
the disk rotation traced in Ha. The molecular gas distribution in NGC 3921 is asymmetric about the 
stellar remnant, and both the distribution and kinematics suggest that the molecular gas has not settled 
into the center of the remnant. Molecular gas clouds are detected outside the central regions of NGC 3921 
and NGC 4676, and they may be associated with the tidal tails and bridges mapped in Hi. Departures 
from the canonical scenario for a merger involving two large spiral galaxies are found in all three Toomre 
Sequence mergers studied. Our data suggest that one of the progentor disks in NGC 520 and NGC 3921 
were relatively gas-poor. A detailed comparison of these optically selected mergers and more luminous 
IR selected mergers is deferred to a companion paper (Paper II). 

Subject headings: galaxies: individual (NGC 520, NGC 3921, NGC 4676) — galaxies: interactions — 
galaxies: evolution — ISM: molecules 



1. INTRODUCTION 

It is widely believed that a physical collision between 
a pair of gas-rich galaxies leads to the concentration of 
gas and an inten se starburst in the central region of the 
merger rern nant (Toomre & Toomre 1972, Negroponte & 



^ ^ [White 1983| ). Indeed the majority of infrared bright galax- 



ies are strongly disturbed systems, and massive concentra- 
tions of molecular gas have been detected in their central 



L'cgions (sec review by Sanders fc Mirabel 1996 and rcfcr- 



oncoD therein). However, thio ia not a complete and gon — 
eral picture of all galaxy collisions since many optically 
disturbed systems do no t show stronglv increas ed levels of 
massive star formation ( Bushouse et al. 1988 ). The ma- 
jority of recent studies have focusscd mainly on infrared 
(IR) selected galaxies, and information on the molecular 
gas content and distribution within merging systems of 
modest infrared luminosity is sparse in comparison. 

A wide range of explanations are possible for the rela- 
tively modest levels of star forming activity observed in 
the optically selected colliding galaxies. Unlike the IR lu- 
minous systems, the progenitor disks of these systems may 
have been relatively gas-poor. Alternatively, the progeni- 
tors may have been gas-rich, but the initial conditions of 
the collisions were such that the bulk of the gas may have 
turned into stars or disp ersed during the collision (e.g. 
Mihos & Hernquist 1996 ). It is also possible that the IR 



As a first step in investigating these possible scenarios, 
we have used Owens Valley Radio Observatory (OVRO) 
Millimeter Array to map the 2.6 mm CO J = 1 tran- 
sition line emission in three "Toomre Sequence" mergers: 
NGC 520, NGC 3921, and NGC 4676. The "Toomre Se- 
quence" is an optically selected ensemble of strongly in- 
teracting galaxies representing a suggested evolutionary 
sequence of disk-disk mergers, based on their stel lar tidal 
tail le ngths and the separation of the two nuclei ( Toomre 
1977). Since the members of this sequence were selected 



luminous phase is generic but brief so that the less lumi- 
nous mergers are seen in the pre- or post-burst phase in 
their evolution. 



on the basis of optical morphology alone, they are much 
less biased towards systems with very high star formation 
rates than IR selected samples. 

The new high resolution interferometric CO observa- 
tions allow us to map the distribution and kinematics of 
the molecular gas in order to investigate the response of the 
molecular material which was previously distributed in the 
inner disks of the progenitors. These data are compared 
with existing high resolution interferometric CO observa- 
tions of IR luminous mergers in a separate paper (Yun 
& Hibbard 2000, hereafter Paper II). In particular, we 
address whether any systematic difference exists in the 
properties of the molecular gas (which directly fuels the 
starburst activity) between the less luminous mergers and 
the IR selected mergers that may explain the differences 
in the IR (and total) luminosity. 

In §2 we describe the three Toomre Sequence mergers 
selected for this study. The details of the observations 
and data reduction are discussed in §3. The results are 
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Table 1 
Summary of Observations 



NGC 520N 



NGC 520S 



NGC 3921 



NGC 4676 



Phase Center 

a (1950) 

5 (1950) 
Center Velocity (km s^"^) 
Gain Calibrator 
Flux Calibrator 
Observed Dates 



On-source Integration (hrs) 
Sensitivity (mjy beam^^) 
9 (FWHM) (PA) 



01:21:58.8 
+03:32:17.0 
2260 
0106+013 



01:21:59.4 
+03:32:09.0 

2260 
0106+013 



11:48:28.9 
+55:21:10.0 
5838 
1150+497 



11.8 
12.9° 



9/28/94, 10/3/94 
12/6/94, 1/26/95 
2/4/95 



2'.'4 X 2''7 (-83°) 



Uranus, Neptune, 3C 273, 3C 454.3 
10/21/94 
2/13/95 
2/20/95 
11.8 19.2 
12.9°^ 5.8*= 

2''3 X 2'.'7 (83°) 



12:43:44.5 
+31:00:02.0 
6600 
1308+326 

10/11/94 
12/5/94 
1/12/95 
16.3 
8.1*' 

2''8 X 3''2 (87°) 



°as measured in each 8 MHz (20.8 km s"^) channel maps, 
^as measured in each 16 MHz (41.6 km s^^) channel maps. 



described individually for each system in §4, and a dis- 
cussion of these results follows in §5. Finally, in §6 we 
summarize our conclusions. 



2. SAMPLE SELECTION 



The Toomre Sequence consists of "Eleven NGC 
Prospects for Ongoing Mergers" , as sketched and pre- 
sented by Toomre (1977). The sequence was arranged 
based on the results of simple numerical simulations con- 
ducted by Toomre & Toomre (1972), which demonstrated 
the effects of bound gravitational interactions on the outer 
regions of disk galaxies. This seminal work illustrated that 
long filamentary features are the natural consequence of 
the tidal forces experienced during such encounters, at 
the rate of one tail per prograde disk. Although their 
numerical technique did not allow for the inclusion of or- 
bital decay, Toomre & Toomre (1972) posited that such 
decay should occur, leading to the eventual merging of 
the participants. This proposal has been repeatedly con- 
firmed with more sophistic ated n umerical treat ments (see 
Barnes & Hernquist 1992, 1996; Barnes 1998 and refer- 
ences therein). The natural consequence of such encoun- 
ters, the Toomres hypothesized, would be pairs of galaxies 
spiraling ever closer together, eventually leaving a single 
stellar body with two protruding tidal tails. The Toomre 
Sequence was meant to depict this proposed evolutionary 
sequence of disk-disk mergers. 

Four members of the Toomre Sequence have previ- 
ously been mapp ed in CO with the OVRO interferometer 
(NGC 4038/ 9 by Stanford et al.^^ , [Wilson et al. 2000|, 
NGC 520 by JSandcrs et al. 1988|; NGC 2623 by [Bryant & 



ScoviUe 1999| ; and NGC 7252 by |Wang et al. 1992| ). These 
systems span the entire range of the sequence from begin- 
ning to end. We chose to improve our understanding of 
the dynamical effects occurring during merging encounters 
by observing additional systems at the beginning and end 
of the sequence (NGC 4676 and NGC 3921, respectively) 
and by re-observing NGC 520. We selected these three sys- 
tems because there are VLA H i mapping and deep broad- 
and narrow-band optical observations available to help in 



the interpretation of the CO data (from Hibbard & van 
Gorkom 1996; hereafter HvG96). We chose to re-observe 
NGC 520 to target explicitly the second nucl eus, which did 



not fa ll within the previously observed field (Sanders et al 



1988). Previous observations and a detailed description of 



the three systems are found in HvG96. Here we summarize 
their defining optical morphological features and their lo- 
cation along the Toomre Sequence. Excellent photographs 
of the syst ems can be found in the Arp "Atlas of Peculiar 
Galaxies" ( Arp 196"^ ), and CCD images of the entire sys- 
tems, included the extended tidal features, are given in 
HvG96. In this paper we will restrict our figures to the 
inner few kpc of each merger, as this is where the CO is 
concentrated. 

NGC 4676 (="The Mice"=UGC 7938/9= Arp 242= VV 
224) is the second member of the Toomre Sequence, rep- 
resenting an early stage merger where the two disks are 
separated by less than one optical diameter in projection 
but still distinct. A bright optical bridge connects the two 
disks, and two distinct bright optical tails, each about 50 
kpc in projected length^, are present. 

NGC 520 (=UGC 966=Arp 157=VV231) falls near the 
middle of the Toomre sequence (7th of 11), representing 
an intermediate stage of merging where two distinct nuclei 
are seen embedded within a single stellar body. The pri- 
mary nucleus is hidden beneath the prominent dust lane 
near the remnant body center, and the secondary nucleus 
is seen 40" (6 kpc) away towards the northwest. A bright 
optical tail stretches 25 kpc to the south, and a broad 
stellar plume is seen reaching 60 kpc to the northwest. 

NGC 3921 (=UGC 6823=Arp 224=Mrk 430) is the next 
to last member of the Toomre Sequence. It represents the 
latest stages of a merger, with a single stellar remnant 
body exhibiting an r^/^ radial light profile, characteristic 
of normal ellipticals. The optical isophotes are not con- 
centric (Schweizer 1996; hereafter S96), suggesting that 
the merger is not fully relaxed. A stellar tail stretches 65 
kpc to the south, and a broad optical plume reaches 80 
kpc to the northeast. 



^AU length, mass, and luminosity scales in the following are derived assuming a value of the Hubble constant of Ho= 75 kms ^ Mpc ^. 
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3. OBSERVATIONS 

Aperture synthesis CO observations of NGC 520, 
NGC 3921, and NGC 4676 were carried out with the 
Owens Valley Millimeter Array between September 1994 
and February 1995. There are six 10.4 m diameter tele- 
scopes in the array, providing a field of view of about 1' 
(FWHM) at 115 GHz. The telescopes are equipped with 
SIS receivers cooled to 4 K, and typical single sideband 
system temperatures were between 300 and 500 K. Base- 
Hnes of 15-200 m E-W and 15-220 m N-S were used, and 
the details of the observations including the synthesized 
beams from naturally weighted data are summarized in 
Table 1. A digital correlator configured with 120 x 4 MHz 
channels (11.2 km s~^) covered a total velocity range of 
1340 km s^^. Nearby quasars (see Table 1) were observed 
at 25 minute intervals to track the phase and short term in- 
strument gain. Uranus {Tb = 120 K), Neptune {Tb = 115 
K), 3C 273, and 3C 454.3 were observed for the absolute 
flux calibration. The data were calibrated using the stan- 
dard Owens Valley array program MMA ( |Scoville ct al 



1992) and were ma pped and analyzed us ing the imaging 
program DIFMAP ( [Shepherd et al. 1994| ) and the NRAO 
AIPS software system. The uncertainty in absolute flux 
calibration is about 15%, mainly due to the uncertainty in 
transferring the calibration between the sources and the 
flux calibrators.'^ The positional accuracy of the resulting 
maps is better than ~ 0'.'5. 



A detailed study of all three systems has been conducted 
in the optical and 21 cm Hi emission by HvG96: radio 
synthesis observations of 21 cm Hi emission to obtain in- 
formation on the distribution and kinematics of extended 
cold atomic gas; deep broad band R images to delineate 
the underlying stellar distribution and any faint optical 
tidal extensions; and narrow band Ha-|-[N il] observations 
to reveal regions of current star formation. These data are 
compared to the CO distribution and kinematics in the 
present work, and are fully described in HvG96. 

Additional optical and near infrared (NIR) data on 
NGC 3921 were needed in order to examine the dynam- 
ical state of this puzzling merger remnant. Broadband 
BVR observations of NGC 3921 were obtained in Jan- 
uary of 1997 with the University of Hawaii 88" (UH88") 
telescope. The f/10 re-imaging optics were used with a 
Tek2048 CCD, resulting in a plate scale of 0.22" pixer^ 
and a field of view of T.b. The seeing was ~1"2, and to- 
tal exposure times of (1200s, 900s, 600s) were obtained in 
{B,V,R), respectively. The data were calibrated via ob- 
servations of Landolt (1983, 1992) standards observed on 
the same nights, with resulting zero-point errors (Icr) of 
(0.01 mag, 0.02 mag, 0.03 mag) in {B,V,R). 

NIR observations were made at X'-band (A = 2.15/^m, 
hereafter referred to simply as K, Wainscoat & Cowie 
1992), obtained with the QUIRC 1024 x 1024 detector 



^Relative accuracy among the measurements presented here is significantly better, but the absolute uncertainty is offered because it is more 
relevant when comparing these measurements with other measurements. 



Figure 1. CO (1—0) spectra of the Toomre Sequence mergers imaged. They are obtained by summing all detected fluxes in the individual channel 
maps (see Figures 3, 7, & 11). Single Dish spectra from the IRAM 30-m telescope {9 pw h M —'^'^ ) are shown in dotted lines (NGC 520 [Solomon et al. 
1992], NGC 3921 [Combes, unpublished], NGC 4676 [Casoli ct al. 1991]). The shape and the flux levels agree well between our measurements and the 
IRAM 30-m measurements, except in NGC 520 where the IRAM 30-m spectrum is significantly lower, possibly due to a pointing or calibration problem. 
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Figure 2. (a)-(c) Velocity integrated CO (1-0) map of NGC 4676 at 2''8 X 3''2 (1.2 X 1.4 Icpc, PA=87°) resolution superposed over the R-band 
image (left), Ha (middle), and 21 cm Hi (right) for comparison. An angular separation of 20" corresponds to 8.5 kpc at the assumed distance of 88 
Mpe. The FOV of the OVRO dishes at 115 GHz (l' FWHM) is indicated by a dashed circle. The CO contours (white) are drawn at levels of 50, 100, 
200, 400, 800 Mq pe~^, and the Hi contours correspond to column densities of 5, 10, 20, 40, 80 times 1 X 10^^ cm"^. 



on the UH88" telescope in January of 1995. The f/10 
re-imaging optics were used, resulting in a plate scale of 
0.187" pixel"^ and a field of view of 3.'2. The observa- 
tions consist of three 120 sec target-sky pairs, with the 
CCD Idithered by 1.'5 between on-source positions. The 



NIR data are uncalibrated. 

4. RESULTS 

The observed and derived properties of the sample are 
summarized in Table 2 along with other properties of inter- 
est. The far-infrared luminosities^, LpiR, of these three 
merging systems range between (1 — 5) x IQ-'^^Lq, with 
a L fir/ Lb ratio of unity for NGC 4676 and NGC 520 
and about 0.2 for NGC 3921. Therefore these are modest 
starburst systems at best. The molecular gas masses^ de- 
tected in CO range between (1 — 5) x 1O^M0, accounting 
for 40% to 100% of the gas masses inferred from single dish 
measurements. This is within a factor of two of the total 
molecular gas associa ted with our Galaxy (3.5 x IO^Mq; 
Sanders et al. 1984). While the bulk of the molecular 



gas in our Galaxy is located within an annulus of 4-6 kpc 
radius, CO emission in these optically selected merger sys- 
tems is concentrated to the central 2 kpc radius, except for 
NGC 4676B, whose total molecular gas content and dis- 
tribution is not well determined by our data (see below). 
The CO (1-0) spectra of each system are produced by 



summing the detected emission from each narrow band 
channel map and are shown in Figure 1. All of the CO 
(1-0) spectra show line widths and shapes co mparable to 



the s i ngle dish spectra (sho wn in dotted lines; Casoli et al 
1991, ^olomon ct al. 1992| , Combes personal communica- 



tion) . We do not recover the full line flux measured in the 
single dish observations in some cases (see Table 2). Since 
the interferometer lacks spacings shorter than about 10 
meters, one possibility is that there is some extended CO 
emission [6 > 45") that is resolved out in our observations. 
Another possibility is that some of the line emission is lost 
to the limiting surface brightness sensitivity of the observa- 
tions. The peak CO line brightness temperature observed 
is only about 0.3-0.6 K for NGC 3921 and NGC 4676 and 
5.2 K for NGC 520. Since the intrinsic CO line brightness 
should be at least 10-20 K and may be as high as 30-50 
K in starburst regions, the beam filling factor for the CO 
emitting regions must be quite small, less than 10% and 
significantly so in some cases. 

The molecular gas distribution and kinematics in the 
individual systems are discussed in detail below. For all 
three observed mergers, we first compare the CO distribu- 
tion with the R-band, narrow band Hq;+[Nii], and VLA 
Hi distributions from HvG96. We then examine the full 
CO 3-dimensional kinematic information as traced by the 
individual channel maps. From these plots it is clear that 



'^Lpin represents the infrared luminosity in the 40-120 fim band, calculated from the 60 /jm and 100/jm IRAS fl uxes (see Hclou et al. 1988). 
Standard Galactic conversion of Nh^/Ico = 3 X 10^0 cm-2 (Kkms-l) ^ is used (see [Young fc Scoville 199l[). 
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Figure 3. The individual cliannel maps of CO (1-0) emission in NGC 4676 with 2'.'8 X 3"2 (1.2 x 1.4 kpc, PA=87°) resolution are shown in 41.6 
km velocity increments upon a greyscale representation of the Ha image. The brightest peak (79.2 m.Ty beam~^) corresponds to a brightness 

temperature of 0.60 K above the microwave background. The contour levels are —2, +2, +3, +4, +5, +6, +7, & +8 times 9 mjy beam~^ (Ic)- The 
dominant CO source is NGC 4676A, and some CO clumps are also seen in NGC 4676B. Several spatially and kinematically coherent CO features in the 
region bridging the two galaxies (V = 6516—6641 km s~^) suggest presence of molecular gas in this overlap region. 



the CO emitting molecular complexes are well resolved 
spatially and kinematically by the aperture synthesis ob- 
servations presented here. We explore the CO kinematics 
by comparing the mean velocity and velocity dispersion 
maps. The mean velocity field gives an idea of the large 
scale orbital motions, and the intensity-weighted velocity 
dispersion helps identify the local centers of the gravita- 
tional potential or the sites of large peculiar velocities. Fi- 
nally, we plot the kinematic profiles along the major axis, 
and compare these to the 21cm line data and any available 
optical kinematics. 



4.1. NGC 4676 
4.1.1. Molecular Gas Distribution 

We mapped the CO (1 0) emission in NGC 4676 (D = 
88 Mpc) at 2'.'8 x 3'.'2 (1.2 x 1.4 kpc) resolution. The 
array was pointed at the region between the two nuclei 
such that the primary beam (field-of-view) of the array 
includes the main bodies of both galaxies, but very little 
of the tidal tails (see Fig. 2a). Both NGC 4676 A (north) 
and NGC 4676B (south) are detected in emission, and the 
integrated CO emission is contoured in Figure 2. About 
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Figure 4. (a) Velocity integrated CO (1-0) map of NGC 4676 in linear eontours superposed over the Ha image. The eontours are 10, 20, 30, 40, 50, 
60, 70, 80, and 90% of the peak which is 18.4 Jy km {Nh2 ^ 0.1 x 10^^ cm"^). Only weak CO emission is detected along the bar in NGC 4676B. 
The spatial correspondence between the CO and Ha is poor, and extinction by dust associated with the molecular gas (peak Ay ~ 120) offers a natural 
explanation in the edge-on system NGC 4676A. (b) Mean CO emission velocity plotted in contours (in km s~^) superposed over the velocity dispersion 
("second moment") map in greyscale. The linear greyscale range between (white) to 160 km s~^ (black). The velocity gradient is increasing to the 
north in both galaxies, and thus the spin orientation of the collision is in a prograde sense. 



80% of the CO emission is associated with NGC 4676A 
while only weak CO emission is detected along the bar in 
NGC 4676B. The total detected CO flux in NGC 4676A 
is 54 ± 8 Jy km s~^ (1 00% of the flux me asured at the 
IRAM 30-m telescope by Casoli et al. 1991 ) , which corre- 
sponds to a total molecular gas mass of 5.5 x 10^ Mq or 
about twice as much as in our Galaxy. In NGC 4676B we 
recover only 20% of the total single dish CO flux reported 
by Casoli et al. Our brightness sensitivity (ATb — 0.08 
K) and a low beam filling factor may explain at least part 
of the "missing" flux. When the visibility data is tapered 
to ^ 5" resolution, the recovered flux increased to 14 ± 2 
Jy km s^^, or about 40% of the total single dish flux. The 
CO morphology does not change substantially in the low 
resolution maps, however. The undetected single dish flux 
is associated with a distinct spiky spectral feature occur- 
ring between the velocity of 6400-6800 km s^^ in Figure 1. 
This feature may arise from the molecular gas associated 
with the bridging region, poorly represented in our data 
due to a low beam filling factor. Alternatively, this fea- 
ture may be CO emission from NGC 4676A picked up by 
the sidelobe response of IRAM 30-m telescope rather than 
being intrinsic to NGC 4676B. 

The CO emission in NGC 4676A is clearly confined to a 
nearly edge-on disk-like structure with a deconvolved size 
of 1.8 kpc in radius and a thickness (FWHM) of 250 pc 
(see Table 2). A bridge of stars and gas connecting the 
two galaxies is clearly seen in the optical and H i emission 



(see Fig. 2c), and several CO clumps are also found in 
the bridging region (i.e. at velocities between 6516-6641 
km s~^ in Fig. 3) albeit with relatively low S/N. Only 
one such CO clump appears in the velocity integrated CO 
map (Figs. 2a & 4a), just southeast of the main body of 
NGC 4676A, because only the high signal-to-noise ratio 
(> 5(t) features are included in these maps. The presence 
of molecular gas in the bridging region suggests that the 
disruption of the inner disks (i? < 10 kpc) has begun in 
NGC 4676. Most of the CO clumps seen in the channel 
map are unresolved (^ 1 kpc in diameter) with molecular 
gas masses of ~ IO^Mq. They are somewhat larger and 
more massive than the giant molecular clouds (GMCs) in 
our Galaxy, and they may be responding ballistically to 
the tidal disruption. Alternatively, these CO clumps may 
be molecular gas condensations forming within the gaseous 
bridge traced in H i and may represent possible sites of fu- 
ture star formation. 

One striking feature in the distribution of CO emission 
in NGC 4676 is the lack of any correspondence between the 
brighter Ha peaks (presumably tracing the present sites of 
star formation) and the CO peaks (dense gas concentra- 
tions). Even a marginal anti-correlation is seen, particu- 
larly in Figure 4a which displays a wider intensity scale 
for the Ha emission than shown in Fig. 2b. The bright- 
est Ha emission in NGC 4676A is located just outside the 
southern tip of the CO emitting region, and significant op- 
tical extinction in this nearly edge-on disk offers a natural 
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Figure 5. Position- velocity plot of CO emission in NGC 4676A along tlie kinematic (and morphological) major axis is shown in greyscale. The 
molecular gas traces the kinematics of the inner disk with a rising rotation curve with a maximum rotation velocity of 270 km s~^, which flattens at 
a radius of 960 pc. The contours represent the corresponding P-V plot of the 21cm H I emission (from HvG), which is a good indicator of the extent 
of the outer disk with a fiat rotation curve. They both suggest that the kinematics of the inner gas disk (i? < 4 kpc) is relatively undisturbed. The 
S-shaped morphology of the CO emission indicates that the molecular gas does not uniformly fill the disk outside 1 kpc radius and may be confined to 
a pair of tightly wound spiral arms. 



explanation ~ the peak integrated CO flux of 18 ± 3 Jy 
km beam"-'^ corresponds to — 6.1 x 10^^ cm^^ 
and a mean visual extinction of Ay ~ 120 averaged over 
the 3" (1.3 kpc) beam. Even the lowest CO contour in 
Fig. 4a corresponds to Ay ~ 12. 

For the more face-on galaxy NGC 4676B, CO and Ha 
emission do not correspond well either. The brightest Ha 
peak near the center of the galaxy has no associated CO 
emission. This may reflect the large extinction associated 
with the CO emitting clouds. Alternatively, this may in- 
dicate that the life time for young stars may be much 
longer than the cloud dispersion time scale. The observed 
"twin-peak" CO morphology is often seen in other barred 
galaxies, probably due t o a bar-driven dyna mical reso- 
nance within these disks (Kcnney et al. 1992). Figure 2b 
shows that both the CO and Ha bars are offset in the same 
manner from the underlying optical bar (HvG96). Similar 
offse ts are seen in hydrodynam ical simulations of merg- 
ers (Barnes fc Hcrnquist 1991) and provide a means of 



NGC 4676B. The apparent contrast in the CO luminosity 
in the single dish measurements by Casoli et al. (1991) is 
about a factor of two, and some of the CO emission as- 
sociated with NGC 4676B may extend beyond the inner 
disk region mapped by us (see above). One explanation for 
the apparent difference in the gas content is that the pro- 
genitor disk for NGC 4676B had less molecular gas. The 
two merging galaxies have similar total H i content (about 
3 X IO^Mq each; HvG96), but NGC 4676B appears to be 
an earlier Hubble type (SBO/a). In a survey of molecular 
gas content among SO and Sa galaxies, Thronson et al. 
(1989) conclude that typical fractional gas masses in SO's 
and Sa's are about an order of magnitude less than those 
for Sb or Sc spirals, and Young & Knezek (1898) report 
the largest M{H2)/M{HI) ratios among the SO/Sa Hub- 
ble types.® Alternatively, the two progenitor disks started 
off with similar amounts of molecular gas but evolved dif- 
ferently under the tidal disruption because of different in- 



transl srring angular momenturn from the gas to the stars, 
allowing the gas to settle even deeper into the potential. 

Another notable aspect of the CO emission in NGC 4676 
is the apparent contrast in CO luminosity and distribu- 
tion between the two merging galaxies. We detect about 
4 times more CO emission in NGC 4676A with a much 
more centrally concentrated distribution compared with 



ternal structure or diff erent spin-orbit alignment (see Mi 



hos & Hcrnquist 1996). The difference between these two 
scenarios has important consequences for understanding 
how a gas disk responds to a tidal disruption. Our obser- 
vations alone cannot distinguish the two however, and this 
issue needs to be addressed by future numerical studies. 



^However, the scatter associated with the individual Hubble type is substantial in both studies. 
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Figure 6. (a)-(c) Velocity integrated CO (1-0) map of NGC 520 at 2''4 X 2'.'7 (0.34 X 0.41 kpc, PA=-83°) resolution superposed over the R-band 
image (left), Ha (middle), and 21 cm Hi (right) for comparison. The location and FOV of the two pointings are indicated by the dashed curves. The 
CO emission is not aligned with the large scale structure of the optical galaxy, which may be the result of accretion of material with angular momentum 
not aligned with the stellar disk. The CO contours (white) are drawn at levels of 50, 100, 200, 400, 800, 1600, 3200 Mq pc~^. The Hi contours 
correspond to column densities of 5, 10, 20, 40, 80, 160, 320 times 1 X 10"^^ cm~^, with dotted contours indicating Hi absorption. 



4.1.2. Molecular Gas Kinematics 

NGC 4676A shows a distinct and relatively intact molec- 
ular gas distribution with relatively undisturbed kinemat- 
ics. The intensity weighted mean CO velocity is shown in 
contours superposed over the greyscale velocity dispersion 
("second moment") map in Fig. 4b. The velocity disper- 
sion shown in Fig. 4b exhibit a distinct peak at the center 
of NGC 4676A. This arises from the rapid rise in the CO 
rotation curve, and therefore marks the dynamical center 
of this edge-on disk. The major axis position-velocity plot 
(Fig. 5) shows more clearly that the CO emitting molecu- 
lar gas in NGC 4676A is in rotation about its center, and 
an apparent flattening of the rotation curve is hinted by 
the sudden drop in the velocity gradient outside the central 
4"-5" (2 kpc) region. The line-of-sight velocity for the CO 
emission sharply decreases at 2'.'5 (1.1 kpc) radius, form- 
ing an S-shaped feature in the position- velocity plot. As 
shown by the solid lines in Fig. 5, the H i rotation curve re- 
mains flat, and this CO kinematic signature probably does 
not indicate a real drop in rotation velocity and may arise 
from non-circular motions induced by a non-axisymmetric 
potential. A pair of tightly wound molecular spiral arms 
can also exhibit a similar kinematic signature. 

The peak rotation velocities traced in CO, HI, and Ha 
are all about 270 km s~^. The large aspect ratio of the 
CO emitting region suggests that the molecular gas disk 
is viewed within 5-10 degrees of being edge-on, and the 
observed peak rotation velocity should be a good estimate 
for the true disk rotation speed. The derived dynami- 
cal mass*" Mdyn = -^tj^ inside the 1.8 kpc radius is then 

^The dynamical mass estimate depends on the assumed mass distribution, and p simplp rntatinnal giippnrt approximation may result in a 
slight over-estimation if the mass distribution is disk-like rather than halo-like (see Binney & Tremaine 1987). 



2.9 X 10^" Mq, and the molecular gas mass inferred from 
the CO emission (4.9 x 10^ Mq) accounts for about 20% 
of this dynamical mass. This gas mass fraction is on the 
upper end of what is seen in ordinary disk galaxies (a few 
to 25%, see Young & Scoville 1991 and references therein). 

The orbital motion of the molecular gas in NGC 4676B 
is more difficult to determine because of its more face-on 
projection and the patchy CO distribution. The observed 
CO kinematics is consistent with the Ha rotation curve de- 
rived by Mihos et al. (1993) along the bar which shows a 
monotonic velocity gradient consistent with solid body ro- 
tation. Assuming a disk inclination of 45°, the dynamical 
mass inside the 4 kpc radius is about 2 x IO^^A/q. The de- 
rived molecular gas mass fraction is then about 7%, which 
is more typical of undisturbed disk galaxies. 

The systemic velocities of the two galaxies are very simi- 
lar (NGC 4676A slightly more redshifted, by ~ 60 km s""'^), 
suggesting that the two disks are either near their orbital 
apocenter or moving mostly in the plane of the sky. The 
velocity gradient increases to the north in both galaxies, 
consistent with the H i kinematics reported by HvG96 and 
the Ha kinematics measured by Stockton (1974) and Mi- 
hos et al. (1993). The spin vector of the two disk are 
aligned with their orbital motion, and the large tidal tails 
emerging from both disks are n aturally explained by thi s 
prograde spin-orbit resonance (Toomre & Toomre 1972). 
The CO emitting clouds seen in the bridging region be- 
tween the two merging disks have the velocities interme- 
diate between the systemic velocities of the two disks as 
expected. 
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4.2. NGC 520 
4.2.1. Molecular Gas Distribution 

We mapped the CO (1-0) emission in NGC 520 {D = 30 
Mpc) at 2'.'3 X 2'.'7 (0.34 x 0.41 kpc) resolution, which is 
a factor of two improvement over the previous OVRO CO 
map by Sanders et al. (1988; 6" x 5"). Two separate 
fields were obs erved, one centered on ea ch of the near- 
infrared nuclei ( |Stanford fc BalceUs 1990 ), but CO emis- 
sion is detected only around the main nucleus (see Fig. 6a) . 
The 3(7 upper limit for the molecular gas mass within the 
1.3 kpc diameter region surrounding the second nucleus 
is 5 x 1O^M0 assuming a total line width of 60 km s~^. 
This is less than 0.3% of the molecular gas mass associated 
with the main nucleus, and the contrast in the associated 
gas mass is quite dramatic. The absence of an associated 
molecular gas complex is somewhat unusual if this were 
a true stellar nucleus of a late type galaxy, and we will 
discuss this point further below (see §5). 

Nearly all of the 3.0 x 10^ A/0 of molecular gas detected 
in the main nucleus position is concentrated in the 1.0 x 0.4 
kpc (PA=95°) disk mapped by Sanders et al. (1988), co- 
incident with the bright, extended nuclear radio source 



( Condon et al. 199C ). The mean H2 density is about 10"^ 
cm^"^ if the molecular gas is uniformly distributed in a disk 
with 500 pc radius and 100 pc thickness, as if the entire 
molecular complex is a single super-massive giant molec- 
ular cloud. The full resolution data results in a detected 
flux of 285 ± 43 Jy km s~^, which is less than the 325 
Jy km s~^ detected by Sanders et al. (1988). Smoothing 
the data to the 6" resolution of the Sanders et al. recovers 
a total flux of 404 ± 61 Jy km s~^. The lower resolu- 
tion maps do not reveal any additional features in the CO 
maps. Our measured flux corresponding to 58% of the to- 
tal CO flux measured by the 14-m FCRAO telescope in 
its central 45" beam area ( Young et al. 199?^ ) . Solomon et 
al. (1992) report a total integrated line flux of only 211 Jy 
km s~^ within the 22" beam of the IRAM 30-m telescope. 
Therefore it is likely that we recovered most of the flux 
associated with the nuclear molecular disk, but there may 
also be significant systematic errors associated with all of 
these line flux measurements compared. 

The partial mapping of CO emission by Young et al. 
suggests that CO emission extends beyond the nuclear re- 
gion, along the stellar body of the galaxy (PA~ 45°). This 
extended CO emission is not detected, either because it is 



Figure 7. The individual channel maps of CO (1-0) omission in NGC 520 with 2'.'4 X 2'.'9 (0.34 X 0.41 kpc, PA=-83°) resolution are shown in 20.8 
km s~^ velocity increments contoured on the greyscale representation of the Ha image. The cross marks the position of the bright radio continuum 
nucleus. The brightest peak (351 mjy beam^^) corresponds to a brightness temperature of 5.2 K. The contour levels are —3, —2, +2, +3, +4, +6, 
-t-lG, -j-15, & -t-25 times 15 mJy beam~^. In addition to the clear rotation signature about the radio nucleus at PA— 0°, there is some evidence for gas 
associated with the large scale tilted disk (PA— -1-45°) or perhaps gas associated with the nuclear blow-out (see channels for V— 2153 &: 2195 km s~^). 
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Figure 8. (a) Velocity integrated CO (1-0) map of NGC 520 in contours superposed over tlic Ha image. The contour levels arc 2, 4, 6, 8, 12, 20, 30, 
and 50 Jy km s~^, and the peak integrated flux observed is 63.9 Jy km s^^ {Nh2 — 3.0 X 10^"^ cm~^). Only the central 5x5 kpc region surrounding 
the primary nucleus of NGC 520 is shown since CO is not detected elsewhere. This integrated CO map has a remarkably similar appearance to that 
of 1.4 GHz radio continuum map by Condon et al. (1990), while the Ha emission appears to avoid the CO emission. A larger Ha image (shown in 
Figure 6) suggests that this ionized hydrogen emission is dominated by the bipolar wind blown out of the intense nuclear starburst region, and the 
starburst is completely obscured by the dust associated with the dusty torus mapped in CO (peak Av ~ 600). (b) Mean CO emission velocity plotted 
in contours (in km s~^) superposed over the velocity dispersion ("second moment") map in greyscale. The linear greyscale range between (white) to 
120 km s^^(black). The well organized velocity fleld and centrally peaked velocity dispersion suggests that the molecular gas in NGC 520 is in rotation 
about the primary nucleus. 



resolved out by the interferometer or because it has insuffi- 
cient filling factor to be detected by the brightness limit of 
the synthesized beam (ATb — 0.18 K). Evidence for some 
extended, diffuse gas is seen in the channel maps (Figure 7, 
espec ially channels at 2195 kms~ ^ and 2153 kms~^), and 



starburst is found in o ther well studied nu clear starburst 
systems such Arp 220 ( ^coville et al. 1997 ). The observed 
anti-correlation between CO and Ha not only confirms the 
diminished H a emission from mas s ive young stars by ex- 
tinction (e.g . Bushouse et al. 1988, Kcnnicutt 1989, Cram 
et al. 1998D , but it also casts a serious doubt on star for- 



its po iition angle with respect to the nuclear disk suggests 
that the extended gas may be associated with either the 
large scale H i disk or with the gas entrained in a galactic 
superwind emerging along the minor axis (HvG96). 

The anti-correlation between Ha and CO emission is 
even more dramatic in the nuclear region of NGC 520. 
The clear displacement of the CO emission from the Ha 
distribution shown in Figure 8a suggests that the Ha emis- 
sion associated with the nuclear starburst is completely ob- 
scured by the dust within the nuclear gas disk. The peak 
integrated CO fiux of 63.9 Jy km s~^ bearn"^ corresponds 
to a mean column density of = 3.0 x 10^^ cm~^ av- 

erage d over the 2'.^5 (375 pc) beam, and the inferred large found in the nuclear s tarburst region in M82 (Lo et al 



mation rates inferred from optical or UV tracers for such 
starburst systems. 

The appearance of the brightest CO emission at the 
most extreme velocities seen in the major axis position- 
velocity plot (Figure 9) as well as the flatness of the CO 
emission along the major axis in the velocity integrated 
map (Fig. 8a) suggests that the CO emission does not rise 
monotonically inwards but has a central hole, i.e. there is 
a nuclear molecular torus rather than a disk (contrast CO 
contours in Fig 8a with those for NGC 4676 A in Fig. 4a). 
This geometry has some similarity to the molec ular torus 



opticf 1 extinction (Ay ^ 600) is entirely consistent with 
the complete obscuration of the nuclear starburst by the 
CO emitting clouds. The observed Ha emission is likely 
dominated by the starburst driven ionized wind escaping 
along the poles and some scattered light from the nuclear 
starburst region. The 1.4 GHz radio continuum map of the 
nucleus of NGC 520 by Condon et al. (1990) is essentially 
identical in appearance and dimension (5" x 2" PA =93°) 
and is coincident with the CO disk shown in Fig. 8a. This 
is further strong evidence that the vigorous starburst ac- 
tivity is indeed associated with the molecular gas complex 
but entirely obscured. Similarly heavy obscuration of the 



1987, Shen & Lo 1995), and such a central hole may be 



a common feature among nuclear starburst systems (e.g. 
Downes fc Solomon 1998[parilh et al. 19981). 



4.2.2. Molecular Gas Kinematics 

The intensity weighted mean CO velocity is shown in 
contours superposed over the greyscale velocity dispersion 
map in Fig. 8b. The CO emitting molecular gas in the 
nuclear region of NGC 520 shows a smooth, monotonic 
velocity gradient of 0.36 km s^^pc^^, which is a clear 
signature of rotation about its center. The regular in- 
terval of the iso-velocity contours suggests a solid body 
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Figure 9. Position-velocity plot of CO emission in NGC 520 along the kinematic (and morphological) major axis is shown in grcyscalc. The constant 
hncar gradient suggests that the molecular disk in NGC 520 lies within the solid rotation part of the inner disk. The dark contours represent the 
corresponding P-V plot of 21cm Hi emission (dotted lines represent Hi absorption - from HvG96), and the eomparison with the CO emission suggest 
that the rotation curve turns over at the location of the outer radius of the molecular disk at a 500 pc radius with a maximum rotation s peed of 200 
km s^^. The brightest CO emission occurs at both extreme velocities, suggesting that the CO emission decreases towards the center (see §4.2.2). 



rotation, which is more evident in the major axis position- 
velocity plot shown in Fig. 9. The 21cm Hi absorption 
(dotted contours) coincides with that of the CO emission, 
and this suggests the presence of neutral atomic gas in- 
termixed with the molecular gas in the nuclear starburst 
region. 

The CO line widths are represented by the greyscale im- 
age in Fig. 8b. The broadening of the line width towards 
the center of the CO complex suggests that the molecular 
gas is well centered on the galaxy nuclear potential. The 
channel maps and the iso-velocity contours suggest that 
the molecular gas disk (or torus) is well resolved and is 
close to but not exactly edge-on. Assuming an intrinsic 
thickness of 100 pc, we infer an inclination of 70° — 75° 
from the aspect ratio of the CO emitting region. The 
peak rotation velocity traced in CO is about 200 km s^^ 
(210 km s~^ correcting for i = 70°). Then the dynam- 



ical n jass inside 500 pc radius is 4.9 x lO'^ Mq, and the Balcells 1990|, HvG96)7 



be smaller by a factor of two or more. The orbital period 
at 500 pc radius is 2 x 10^ years, which is comparable to 
the time scale for the starburst but two orders of magni- 
tude smaller than the dynamical time scale for the merger 
(HvG96). 

One unusual aspect of the nuclear gas torus is that 
its position angle and angular momentum vector are mis- 
aligned with the larger scale stellar structures and outer 
H I disk. A similar misalignment of the nuclear disk is 
seen in other nuclear starburst systems like M82 and may 
be produced by accretion of inflowing gas with misaligned 
angular momentum. The observed behavior is probably 
transient in nature as the continuing infall of mass and 
resulting torque will continue to shape the potential. The 
systemic velocity of the molecular gas disk, 2247 km s~^, 
is in good agreeme nt with to that of the lar g e scale stel- 
lar and Hi d isks ( ptockton fc Bcrtola 1980 , Stanford & 



H2 mass inferred from the CO emission (4.3 x 10'=' Mq, 
see Table 2) nearly entirely accounts for the total mass 
inside the 500 pc radi us reg ion. This is a larger fraction 
than in NGC 4676A (§4.1.2) and is similarly dominant as 
i n the IR luminous mergers such as Mrk 231 and Arp 220 
( [Bryant fc Scoville 1996| , |Scoville et al. 1997| ). The higher 
density and temperature conditions associated with the in- 
tense starburst regions may result in an over- estimation of 
molec ular gas mass if the standard C0-to-H2 convers ion is 
used (Scoville et al. 1997, Downes & Solomon 199S), and 
the actual molecular gas mass traced by CO emission may 



4.3. NGC 3921 

4.3.1. Molecular Gas Distribution 

NGC 3921 (D = 78 Mpc) was detected previously in 
CO (1-0) and CO (2-1) emission with the IRAM 30 m 
telescope by Combes and collaborators (Combes, personal 
communication). We mapped the CO (1-0) emission in 
NGC 3921 at 2'.'3 x 2'.'7 (0.87 x 1.0 kpc) resolution, finding 
the emission confined to at least two separate complexes 
(Figure 10): a main complex of 2 kpc diameter near the 
center of the field, and an unresolved CO emitting cloud 
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Figure 10. (a)-(c) Velocity integrated CO (1-0) map of NGC 3921 at 2''3 X 2''7 (0.87 X 1.0 kpc, PA=83°) resolution superposed over the R-band 
image (loft), Hq (middle), and 21 cm Hi (right) for comparison. An angular separation of 20" corresponds to 7.5 kpc at the assumed distance of 78 
Mpc. The CO emission appears centered on the remnant on these angular scales (15^'), but it is displaced from the brightest optical peak as seen in 
the comparison with Ha and as in Figure 12. The CO features occur near the ridge of bright HI emission, but they do not coincide with the H I peaks. 
The CO contours (white) are drawn at levels of 50, 100, 200, 400 Mq pc~^, and the Hi contours correspond to column densities of 5, 10, 20, 40, 80 



times 1 X 10 cm 



located 7" (2.6 kpc) south of the nucleus. The inferred 
molecular gas masses are 1.6 x 10^ Mq and 2.0 x 10^ 
Mq, respectively. We recover all of the CO flux detected 
at the central position by the 22" beam of the IRAM 30 
m telescope (see Fig. 1 & Table 2), but the CO emission 
may extend beyond the area we mapped (Combes, per- 
sonal communication). The mean H2 density for the main 
molecular gas complex is ^ 150 cm~^ if the assumed line of 
sight thickness is 100 pc. This mean density is about 1/6 
of that in NGC 520 but comparable to the typical density 
in Galactic CMC's. 

The integrated CO emission map and the derived mean 
velocity field are shown in comparison to the HI and Ha 
emission in Figures 10 & 12. The locations of the optical, 
near-infrared, and radio continuum peak are coincident to 
within a fraction of an arcsecond, and the mean position 
is marked with a cross. The main molecular complex in 
NGC 3921, however, is displaced by about 2" (- 760 pc) 
to the west from this location. This suggests that either 
the gas is physically displaced from the optical nucleus, or 
that the true nucleus lies totally obscured beneath the CO 
complex (as is the case for NGC 520 and NGC 4676A). 

These possibilities are further investigated by compar- 
ing the CO distribution with the optical and NIR (2.15;Um) 
images, as shown in Figure 13. The optical and NIR im- 
ages are found to be coincident to within a fraction of a 
pixel (~ 0.2"), while the peak of the CO emission is dis- 
placed by at least 10 pixels (2") to the west. Further, the 
ivT-band image and B — K color map shows no hint of a 
hidden NIR peak beneath the CO complex. While the 
peak extinction inferred from our CO observation is large 
{Ay ^ 70, Ak ^ 7), evidence for a second nucleus would 



be visible in the NIR image if an extended stellar bulge 
is present. It therefore seems that the observed displace- 
ment between the molecular gas and the stellar nucleus is 
real. The CO emission appears to be associated with a 
concentration of dust, as indicated by the dark greyscales 
in the B — K map in Fig. 13. 

A comparison between the gas distribution and the HST 
Planetary Camera (PC) image of NGC 3921 by Schweizer 
et al. (1996; see especially their Figure 4) supports this 
picture. This comparison is shown in Fig. 14, where the 
CO emission is contoured upon the F-band PC image of 
NGC 3921 of Schweizer et al. and upon the same image af- 
ter a best-fit model light distribution has been subtracted. 
In the latter image dust lanes show up in white, and we 
confirm the conclusion reached above that the main CO 
complex coincides with an intricate system of dust lanes 
concentrated mainly to the west of the nucleus. Such a 
displacement is unexpected since the highly dissipative 
nature of the cold gas should cause it to settle into the 
gravitational potential within an inner dynamical time of 
^ 10^ years. The non-axisymmetric appearance of the 
dust lanes gives the impression that the gas continues to 
spiral in towards the center, and the gas may be still set- 
tling within the merger potential. The co mbination of off - 
centered or "sloshing" optical isophotes ( 3chwcizer 1996 ) 
and non-axisymmetric dust lanes and disordered kinemat- 
ics led Schweizer et al. (1996) to suggest that "on scales 
of order 100-1000 pc the gas and dust in NGC 3921 do 
not form a well settled nuclear disk" . The displaced CO 
appears to be a strong confirmation of this interpretation. 

The molecular gas clump located to the south of the 
main complex lies along the ridge of optical. Ha, and H i 
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Figure 11. The individual ciianncl maps of CO (1-0) emission in NGC 3921 with 2'.'3 X 2'.'7 (0.87 X 1.0 kpe, PA=83°) resolution are shown in 
41.6 km s~^ veloeity increments contoured on a greyscale representation of the Ha image. The brightest peak (42.1 mjy beam~^) corresponds to a 
briglitness temperature of 0.64 K averaged over the beam area. Tlie contour levels are —3, —2, +2, +3, +4, & +5 times 8 mJy beam~^ (1*^). The CO 
emission occurs clearly displaced to the northwest of the radio nucleus, and the gas kinematics is more complex than a simple rotation. 



emission that delineates the beginning of the southern tail 
(see Fig. 10, and also HvG96). This raises an intriguing 
possibility that this feature is associated with the gas-rich 
southern tail. Tidal features are generally very difficult 



to detect in CO (e.g. Smith fc Higdon 1994, Smith et al 



1999^] but molecular gas h as be en found in the bridging 



regio ns in N GC 4676 (see § [l.l.l| ) and in Arp 105 (Brainc 
et al. I 200C ). The continued infall of gas from the outer 
regions might offer an explanation for the displaced molec- 
ular complex in NGC 3921 (see §5). 

4.3.2. Molecular Gas Kinematics 

The CO channel maps for NGC 3921 are given in Fig. 11. 
These maps show a general southwest-to-northeast veloc- 
ity gradient along the main CO complex. However, numer- 
ous peaks and extensions appear and disappear in different 
channels at various position angles. The southwest com- 



plex appears at velocities between 5752-5838 km s^^, but 
with a velocity gradient opposite in sense to the main com- 
plex (i.e., decreasing velocities from south to north). The 
intensity weighted mean CO velocity, shown as contours 
superposed over a greyscale of the velocity dispersion map 
in Fig. 12b, also show the same gradient along the main 
CO complex and the velocity reversal in the SW complex. 
On the other hand, the local maximum in velocity disper- 
sion characteristic of a local gravitational potential well is 
not seen in Fig. 12b, supporting the conclusion reached 
earlier that no distinct potential maximum exists within 
the main molecular gas complex. 

Towards the main body, the H i column density drops 
rapidly, and the H i tail cannot be traced as an individ- 
ual kinematic structure. The H i which appears projected 
onto the main body is spread over a broad range of veloci- 
ties, similar in width to the OVRO and single dish profiles 



14 



Molecular Gas in Optically Selected Mergers 



Table 2 
1 



and the stellar ve locity disper- 
Lake fc Dressier 1986|). Unlike the 



(~ 460 km s 
sion {r-^ 280 km s 
case for the other two systems, the molecular gas mass 
accounts for a relatively small fraction of the dynamical 
mass in the inner regions (about 10% and 5% inside a 1 
and 2 kpc radius, respectively). 

The position-velocity plot of the CO data is not very 
informative, so we instead illustrate the CO kinematics by 
comparing them to the optical emission line kinematics of 
S96, which were taken through a 2" slit centered on the 
optical nucleus at a PA of 45°. This comparison is pre- 
sented in Figure 15. The emission line kinematics are in 
general agreement with that the CO, suggesting that the 
two are tracing the same feature. To the SW, both the 
optical emission lines and CO show a sudden velocity re- 
versal to lower velocities, showing that neither component 
clearly traces a simple rotation about the central poten- 
tial. S96 derives a systemic velocity of 5926 km s~^ from 
the emission lines, while the CO data may suggest a some- 
what lower value. However, given the lack of a recogniz- 
able systematic rotational signature, we believe the data 
are inconclusive as to this point. 

The observed velocity reversals are very similar to those 
seen in the well studied mer ger remnant NG C 7252. A 



comparison between the Ha (Schwcizcr 1982), CO ( Wang 



1992|), and t idal Hi kinematics in th at system^ Hi 



et al. 1994, Hibbard fc Mihos 1995) show that the 



et al 
bbard 

velocity reversals occur where tidal features are projected 
near the inner regions, and can be explained by suppos- 
ing that strong streaming motions of infalling gas from 
the tidal regions are superposed upon the kinematics of 
a central molecular disk. In NGC 3921 the velocities of 



the southern clump (encircled "S" in Fig. 15) appear at 
similar velocities as the blueshifted ionized gas. These ve- 
locities are also similar to the H i velocities at the base of 
the southern tail (HvG96), and this kinematic component 
may be associated with infalling gas from the tidal regions. 

In summary, the inner region of NGC 3921 has an off- 
center molecular gas complex which has not settled into a 
disk, and a second component of gas may be infalling from 
the tidal regions. This is rather different from the nuclear 
molecular complexes seen in most single nucleus mergers, 
a point we will return to below. 

5. DISCUSSION 

The majority of mergers mapped in CO show a cen- 
tral molecular disk with well defined rotational kinematics 



(Scovillc et al. 1997, Downcs fc Solomon 1998, Bryant fc 



Scovillc 1999| ). Early gas inflow and the formation of a 
compact nuclear gas complex, even before the coalescence 
of the two stellar nuclei, are also seen in numerical studies 
of galaxy merger (e.g. Barnes fc Hcrnquist 199(^ ). This 
is a direct consequence of the fact that gas can shock and 
radiate energy away. 

While some of the observed properties of these three 
Toomre Sequence mergers follow this canonical picture of 
a merger involving two large spiral galaxies, significant de- 
partures are also found. A dense concentration of molec- 
ular gas is associated with the central nucleus of the late 
stage merger NGC 520, but no other molecular gas con- 
centration is detected on the second nucleus or along the 
region bridging the two nuclei. The formation of a com- 
pact nuclear gas complex appears to be well underway in 
the merging disk of NGC 4676A, but evidence for gas in- 



Figure 12. (a) Velocity integrated CO (1-0) map of NGC 3921 in linear eontours superposed over the Ha image. The contours are 10, 20, 30, 40, 
50, 60, 70, 80, and 90% of the peak which is 6.4 Jy km s"^ (JVffs = 3.7 X 10^^ cm~^). The bulk of the 1.8 X 10^ Mq molecular gas is concentrated 
in a single 2 kpc diameter complex, which is clearly displaced from the optical and radio continuum peak (marked with a cross). The HST images by 
Schweizer et al. (1996) reveal an abrupt change in color and an extensive array of dust filaments in the same region where the CO emission is detected 
(peak Av ~ 70 — see Fig. 13). The displacement of molecular gas from the center of the optical remnant appears real and not an extinction effect, (b) 
Mean CO velocity plotted in contours (in km s^^) superposed over the velocity dispersion ("second moment") map in greyscale. The linear greyscale 
range between (white) to 150 km s~^ (black). Unlike NGC 520 or NGC 4676A, there is no clear potential center characterized by a large velocity 
dispersion (AV > 100 km s^^), and evidence for reversals in the velocity field is seen. 
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Figure 13. Comparisons of integrated CO emission in NGC 3921 with B, K , and B — K images. A cross in eaeh image marks tlie center of NGC 3921 
as determined from tlie HST observations of Schweizer et al. (1996). The spatial offset of the central CO complex from the stellar nucleus is apparent 
in all continuum band images. Absence of any K-band emission coincident with the CO feature rules out the possibility of an obscured nucleus within 
the molecular gas complex. A spatial correspondence is found between the CO emission and the dust lanes to the northwest (darker B — K colors in 
the lower panels). The CO is contoured at 10, 20, 90% of its peak value. North is up and East to the left, with vertical tic-marks drawn every l" 
and horizontal tic-marks drawn every 0.1 sec. 



flow and formation of a central gas concentration is not 
evident in NGC 4676B. In NGC 3921, the molecular gas 
complex is significantly displaced from the peak of the op- 
tical, near-infrared, and radio continuum emission 760 
pc). 

One possible explanation for these departures from the 
canonical merger scenario is that one or both progenitor 
disks of these optically selected merger systems may have 
been relatively gas-poor. Based on the analysis of the stel- 
lar and H I tidal features, Hibbard & van Gorkom (1996) 
suggest that both NGC 520 and NGC 3921 are the prod- 
ucts of an encounter between one gas-rich disk and one 
gas-poor system, such as an SO or Sa galaxy (e.g. Thron- 
son et al. 1989). 

This suggestion for NGC 3921 is further supported by 
the observations of Schweizer et al. (1996), who find many 
fewer young globular clusters in this system compared to 
the gas rich mergers NGC 4038/9 and NGC 7252. The ab- 
sence of molecular gas directly associated with the second 
nucleus of NGC 520 and the stellar nucleus of NGC 3921 



is then naturally explained if little molecular gas has been 
tunneled into the nuclei during the merger process. The 
absence of gas in one progenitor also lessens the frequency 
of gas cloud collision s in the second disk, e specially within 
the central few kpc ( Olson fc Kwan 1990 ), and therefore 
less dissipation and angular momentum transport is ex- 
pected. 

The presence of two large, atomic gas-rich tidal tails 
associated with NGC 4676 suggests that both progenitor 
disks were gas-rich (see HvG96), and a different expla- 
nation may be needed. As stated earlier, NGC 4676B 
appears to be an early Hubble type spiral and may also 
have been relatively molecular gas poor initially. A nu- 
merical study by Mihos & Hernquist (1996) has shown 
that the presence or absence of a massive bulge in the 
progenitor disk can produce a large difference in the evo- 
lution of the gas distribution and kinematics, the peak gas 
density achieved, and the resulting star formation activ- 
ity. The apparent difference between the current molecu- 
lar gas content in the two merging disks in NGC 4676 may 
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Figure 14. Comparisons of integrated CO emission in NGC 3921 witli (left) tlie V-band HST image from Scliwcizer et al. (1996, tlieir Figure 4a), 
and (right) witii tiie same image, after a best-fit model liglit distribution has been subtracted (their Figure 4b). The region shown is approximately the 
same as in Fig. 13, although the field-of-view has been rotated. The model-subtracted image shows much more clearly the dust lanes concentrated to 
the west of the system, which appear white in this figure. The CO is contoured as in Fig. 13, and is seen to be very closely associated with the dust 
lanes to the west. The chaotic appearance of the dust lanes led Schweizer et al. (1996) to suggest that NGC 3921 is still in the process of forming a 
nuclear disk. The displaced CO appears to be a strong confirmation of this interpretation. 



be explained more naturally by differences in the internal 
structure of the progenitors. 

While it is significant that these three optically se- 
lected mergers show notable departures from the canoni- 
cal merger scenario, we should also keep in mind that they 
nevertheless represent only three individual snapshots of 
their respective merger evolutions. For example, the ap- 
parent difference in the status of molecular gas between 
the two merging disks in NGC 4676 may simply reflect a 
slight delay in the onset of the gas inflow created by the 
difference in the spin-orbit coupling, however fortuitous 
this may seem. Similarly, NGC 3921 may be in a brief, 
early stage in its merger evolution and may soon develop 
a molecular disk. Schweizer (1996) has also presented ob- 
servation evidence that NGC 3921 is a dynamically young 
merger. However, central CO disks are commonly found 
in both earlier (e.g. NGC 520 and NGC 4676A) and later 
stage (e.g. NGC 7252) mergers. One may further postu- 
late that NGC 3921 represents a stage even further in its 
evolution where the original central gas complex is already 
dispersed. In this scenario, the observed gas in the central 
region is largely the result of recent inflow from the tidal 
tails (see Hibbard & Mihos 1995). 

Finally, a clear and significant quantitative result from 
this imaging study of optically selected merger systems 
is that their central molecular gas surface density is sys- 
tematically smaller by orders of magnitude compared with 
the IR selected mergers (see Table 2). Further, the level 
of activity associated with each of the merger systems ap- 
pears to depend strongly on the central gas density. We 
will discuss this and other comparisons of global properties 
between the optically selected and IR selected mergers in 
a separate paper (Paper II). 



6. SUMMARY 

In an effort to better understand the nature of the lu- 
minous infrared galaxy phenomenon and to scrutinize the 
details of the merger scenario, we have mapped the CO 
emission in an optically selected sample of on-going merg- 
ers with mo derate IR luminosity ( "Toomre sequence" , 
Toomre 1977) at a resolution and sensitivity comparable to 
those of the more IR luminous objects, and the response of 
the molecular gas material deep inside the potential of the 
colliding galaxies is traced. The summary of our findings 
is as follows: 

1. A compact (i? < 2 kpc) molecular complex is found 
well centered on the inner stellar disk of NGC 4676A, form- 
ing a disk- or a ring-like structure with the regular kine- 
matics consistent with simple rotation. This molecular 
gas concentration contributes a significant fraction (20%) 
of the total mass in the nuclear region. The CO emission 
occurs along the 7 kpc stellar bar in NGC 4676B. The sur- 
face brightness of the CO emission is extremely low, and 
tapering the data from 3" to 5" resolution nearly doubles 
the detected flux. 

2. Nearly all of the CO emission detected at the primary 
nucleus position in NGC 520 is concentrated in a 1 kpc di- 
ameter ring-like structure, coincident with the extended 
radio continuum source mapped by Condon et al. (1990). 
The derived H2 mass and the dynamical mass are compa- 
rable, and the gas mass must constitute a large fraction 
of the total mass in the nuclear starburst region. No CO 
emission is detected on the second nucleus in NGC 520. 

3. The molecular gas distribution in NGC 3921 is quite 
different from any other singly nucleated merger mapped 
in CO. The CO emission is significantly displaced from the 
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Figure 15. Position-velocity plot of CO emission in NGC 3921 in comparison to the kinematics of the optical emission lines is shown on the right 
panel while the actual positions of each sampled component are shown on the right panel. The Ha kinematics of S96 arc indicated by filled circles. CO 
velocities measured from the first moment map across the peak of the CO distribution at a PA of 45° are shown as small 's. The "clean components" 
associated with the southern clump arc denoted by encircled "S" . The size of the symbol is proportional to the S/N of the component. 



peak jf the optical, near-infrared, and radio continuum fc Hcrnquist 1996| ) may offer a plausible explanation for 



emission (~ 760 pc), and both the distribution and kine 
matics suggest that the molecular gas has not settled into 
the new remnant potential. The molecular gas complex 
contributes a minor fraction 10%) of the dynamical 
mass in the inner regions. 

4. Molecular gas clouds are detected outside the cen- 
tral regions of NGC 3921 and NGC 4676. They may be 
associated with the tidal tails and bridges mapped in H i. 

5. A consistent trend of anti-correlation is seen between 
CO and Ha emission in NGC 520 and NGC 4676A, and the 
large extinction inferred from the CO emission {Ay = 600 
& 120) suggests that the intense nuclear starburst regions 
in these galaxies are completely obscured. This finding of- 
fers a natural explanation for the severe under-estimation 
of the star fo rmation rate by optica l tracers compared with 
the IR (e.g. [Bushouse et al. 1988| , Kennicutt 1989,1998). 

6. Departures from the canonical scenario for a merger 
involving two large spiral galaxies are found in all three 
Toomre Sequence mergers studied. A relatively gas-poor 
progenitor, as inferred from the optical and H i observa- 
tions by HvG96, offers a possible explanation for NGC 520 
and NGC 3921. A gas-poor progenitor disk or rapid evolu- 
tion in gas content by some dynamical process (i.e. Mihof 



the low molecular gas density in NGC 4676B. 

7. The central molecular gas density in these optically 
selected mergers is systematically smaller by an order of 
magnitude compared with the IR selected mergers, the 
level of activity associated with each of the merger sys- 
tems appears to depend strongly on the central gas den- 
sity. This and other detailed comparisons of global prop- 
erties between the optically and IR selected mergers will 
be presented in Paper II. 
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Table 2 

Summary of Observed and Derived Properties 





NGC 4676A 


NGC 4676B 


NGC 520 


NGC 3921 


RA (1950) 


12:43:44.1 


12:43:45.3 


01:21:59.6 


11:48:28.9 


DEC (1950) 


+31:00:19 


+30:59:47 


+03:31:53 


+55:21:23 


Distance (Mpc)^ 


88 


88 


30 


78 


< Vco.hcl > (km s"^) 


6632±42 


6590±42 


2247±21 


5880±42 


^Vco.FWZl (km s^'-) 


546±42 


252±42 


499±21 


462±42 


Dcconvol. Size (PA) 


8''5 X I'.'l (5°) 


16" X 2"(22°) 


7''1 X 2''9 (95°) 


5''0 X 2''3 (38°) 


(3.6 X 0.5 kpc) 


(7x < 1 kpe) 


(1.0 X 0.4 kpc) 


(1.9 X 0.9 kpc) 




0.60 ± 0.09 K 


0.28 ± 0.04 K 


5.2 ± 0.7 K 


0.64 ± 0.09 K 


Sco^V (Jy km s'^) 










OVRO 


54 ± 8 


14 ± 2 


404 ± 61 


25 ± 4 


Single Disli"" 


55 


34 


696 


21 


Mh^,OVRO (Mq)" 


(4.9 ± 0.7) X 10 


' (1.3 ± 0.2) X 10^ 


(4.3 ± 0.7) X 10^ 


(1.8 ± 0.3) X lO'' 


(6.1 ± 0.9) X 10^^ (1.6 ± 0.0.2) X 10^^ 


(3.0 ± 0.5) X 10^^ 


(3.7 ± 0.6) X 10^^ 


(Mo pc-2)'^ 










peak 


900 ± 135 


210 ± 32 


4670 ± 700 


540 ± 270 


mean 


340 ± 51 


24 ± 4 


2650 ± 400 


134 ± 20 


Si. 4 Gffi (mjy)° 


23.6 


6.6 


158 


9.4 


log Li,4 GH. (W Hz-i) 


22.0 


21.4 


21.9 


21.5 


IRAS 12/i (Jy)*" 




(0.09 ± 0.03) 


0.77 ± 0.06 


0.12 ± 0.02 


IRAS 25/i (Jy)*" 




(0.48 ± 0.03) 


2.87 ± 0.17 


0.28 ± 0.02 


IRAS 60m (Jy)*" 




(2.93 ± 0.04) 


30.9 ± 1.9 


0.83 ± 0.03 


IRAS lOO/i (Jy)*" 




(5.26 ± 0.12) 


45.8 ± 2.7 


1.60 ± 0.11 






38 K 


42 K 


37 K 


Lfir [Lq)^ 


2.8 X 10^" 


0.8 X 10^" 


5.4 X 10^" 


0.8 X 10^° 


Lb {W^°LqY 


2.8 X 10^° 


1.6 X 10^° 


3.4 X 10^" 


5.5 X 10^° 



Ho — 75 km s ^ Mpc ^ assumed. 

Single Dish references: NGC 4676: Casoli et al. (1991); NGC 520: Young et al. (1995); NGC 3921: Combes, personal communication. 

H2 



= 1.2 X lO^'ScoAVD^jp^ Mq (see Sanders et al. 1991) 



M 



TVR'^ 1-1 
ln2 1 



H2 

"From HvG96. 

'from IRAS ADDSCAN/SCANPI (see Holou et al. 1988). 



^Lfir — 4 X 10"^ [2. 58560/^ + S-ioOtA^ Mpc assuming dust cmissivity n — 1 (see Helou et al. 1988). NGC 4676 A&B are not 

resolved by IRAS, and the FIR fluxes are estimated by the ratio of their 1.4 GHz radio continuum flux ratio. 
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